A scalar two-level boson model to study the IBM phase diagram in the Casten triangle 



O 
O 

>>'. 



(N 
> 

\o 

O 

(N 
O 
VO 
O 



o 

G 



X 



Julien Vidal, 1 Jose M. Arias, 2 Jorge Dukelsky, 3 and Jose Enrique Garcfa-Ramos 4 

1 Laboratoire de Physique Theorique de la Matiere Condensee, CNRS UMR 7600, 
Universite Pierre et Mane Curie, 4 Place Jussieu, 75252 Paris Cedex 05, France 
2 Departamento de Fisica Atomica, Molecular y Nuclear, Facultad de Fisica, 
Universidad de Sevilla, Apartado 1065, 41080 Sevilla, Spain 
Institute- de Estructura de la Materia, CSIC, Serrano 123, 28006 Madrid, Spain 
* Departamento de Fisica Aplicada, Universidad de Huelva, 21071 Huelva, Spain 

We introduce a simple two-level boson model with the same energy surface as the Q-consistent 
Interacting Boson Model Hamiltonian. The model can be diagonalized for large number of bosons 
and the results used to check analytical finite-size corrections to the energy gap and the order 
parameter in the critical region. 
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I. INTRODUCTION 

The study of quantum phase transitions is a hot 
topic covering different branches of quantum many-body 
physics. In recent years, the study of low-dimensional lat- 
tice models has been triggered by the increasing interest 
on the development of a quantum computer. Moreover, 
concepts of quantum information have been used to char- 
acterize quantum critical phenomena. On a different re- 
spect, there is a revival of the study of structural changes 
in finite-size systems where precursors of the transitions 
could be observed 0. In atomic nuclei, the Interacting 
Boson Model (IBM) provides a framework simple but 
still detailed in which first and second order phase tran- 
sitions can be studied. In fact, the transition from the 
spherical t/(5) dynamical symmetry to the deformed 7- 
unstable 0(6) dynamical symmetry has been extensively 
studied in the last few years 0, A H| • The Hamilto- 
nian describing this transition is a repulsive boson pairing 
Hamiltonian that has the particularity of being exactly 
solvable allowing the study of very large systems. The 
boson pairing model was first solved by Richardson |6| 
and more recently applied to the IBM However, 
this is the only integrable line in the phase diagram of the 
IBM often described with the Q-Consistent Hamiltonian 
(QCH) ^3 and pictorially represented by the Casten tri- 
angle. This integrable line corresponds to one side of the 
triangle. Therefore, the interior area of the triangle as 
well as the other two sides are out of reach of large scale 
numerical studies. Numerical calculations in these cases 
are then limited by the number of bosons that standard 
IBM codes can manage, which are of the order of 10 2 
bosons. However, strictly speaking quantum phase tran- 
sitions only occur in macroscopic systems. Thus, it is 
of great interest to be able to extend previous analysis 
in the Casten triangle to a significant larger number of 
bosons. 

In this paper, we introduce a simple two-level boson 
model depending on two control parameters, that leads to 
the same energy surface as the Q-consistent IBM Hamil- 
tonian when the 7 degree of freedom is frozen at 7 = 0. In 



other words, both energy surfaces coincide if axial sym- 
metry is imposed in the IBM treatment. This model can 
be solved for very large boson numbers and it allows us to 
investigate the properties of the first-order phase transi- 
tion from spherical to deformed axial shapes in the IBM. 
The two-level boson model can be viewed as a general- 
ized Lipkin model with parity breaking terms, as a two 
coupled large spin system or as a model to study 
tunnelling dynamics between two minima |12|. 

In Sect. II, after a brief presentation of the model, we 
derive its energy surface. Next, in Sect. Ill we compute 
the finite-size corrections up to order (1/iV) 1 in the spher- 
ical phase which allows us to discuss the finite-size be- 
haviour at the critical point for different quantities such 
as the ground state energy, the gap, and the order pa- 
rameter. Finally, Sect. IV is for summarizing. 



II. THE MODEL 

Let us now introduce the two-level boson model whose 
Hamiltonian is 



where the operators n t and Q v are defined as 
n t = th, Q v = sH + tfs+y th 



(1) 



(2) 



in terms of two species of scalar bosons s and t, x and 
y being two independent control parameters. The total 
number of bosons N = n s + n t is a conserved quantity. 
The connection between the Hamiltonian Q and a gener- 
alized Lipkin model or a two-spin model can be obtained 
by making the inverse of the Schwinger transformation 
K+ = tU, K- = sH, and K° = \{th - sU). 

We have deliberatively written the two-level boson 
Hamiltonian in the form to resemble QCH, where 
s and t play the role of the s and bosons of the 
IBM respectively. Clearly, the difference resides on the 
quadrupolar character of the d M boson leading to the 
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U(6) algebra of the IBM, while our model is described by 
a U(2) algebra. Despite of this important difference, the 
st model captures the main characteristics of the phase 
diagram of IBM as we will now show. 

The phase diagram of model (JTJ can be easily obtained 
in the coherent state approach. Therefore, we introduce 
a variational state of the form 

\N,p) =e ^TT^( s>+l3ti ) | ) (3) 

where |0) denotes the boson vacuum. The corresponding 
energy surface as a function of the variational parameter 
(3, in the large N limit, is given by: 

E(N,0) = (N,0\H\N,fl 

f3 2 [x + y 2 (x-l)]}. (4) 

Minimization of the energy Q with respect to (3, for 
given values of the control parameters x and y, gives the 
equilibrium value (3q defining the phase of the system. 
/3 = corresponds to the symmetric phase, and 0o ^ 
to the broken symmetry phase. Here, we restrict the 
study to the range x € [0, 1] and y > which implies 
A) > 0. However, the treatment can be easily extended 
to the whole parameter space. 

The schematic phase diagram of the st model is pre- 
sented in Fig. \I]m the region of interest around the phase 
transition. In addition to the critical line displayed by 
the solid curve, we show the spinodal and antispinodal 
lines where the second minimum in the energy surface ijljl 
starts to appear indicating a region in which both phases 
coexist. 
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FIG. 1: (Color online) Phase diagram for the Hamiltonian Q 
in the plane of the control parameters (x,y). For y = 0, the 
transition is second-order whereas it is first-order otherwise. 

For y — 0, there is an isolated point of second or- 
der phase transition as a function of x. Spinodal, anti- 
spinodal and critical point coincide at the critical value 



x c = 4/5. For y ^ the phase transition changes its 
character to first order as can be seen in Fig. [21 for 
y = 1/ a/2- Schematically, for x = 1, the system is in the 
symmetric phase since the energy surface has a unique 
minimum at f3 — 0. When x decreases, one reaches the 
spinodal point [x = 0.82036 for y = 1/V2) where a sec- 
ond local (nonsymmetric) minimum {(3 ^ 0) arises. This 
nonsymmetric minimum competes with the symmetric 
one till both attain the same depth at the critical point 
x c = (4 + y 2 )/(5 + y 2 ). Beyond this value, the sym- 
metric minimum at f3 = becomes a local minimum till 
x = 4/5 where it becomes unstable (anti-spinodal point). 
We show in Fig. (left panels), a sketch of this evolu- 
tion for the special cases y = 0, 1/V2- The right panels 
in Fig. [5] show two cases in the coexistence region. 

It is worth mentioning that this scheme is exactly the 
same as the one obtained in the IBM of nuclear structure. 
In fact the value y = 1/V2 is equivalent to x = — v7/2 
describing the t/(5) to SU(3) side of the Casten triangle 
[l3| as we are showing now. The energy coincides 
with that of the Q-consistent IBM-1 Hamiltonian in the 
thermodynamic limit (large N). This correspondance is 
readily established if we write the Q-consistent IBM-1 
Hamiltonian 

1 — x 

Hibm = x n d — Q x ■ Q x , (5) 

where rid is the d— boson number operator and Q x is 
defined as 

(2) 

Q x = stj + dKs + x (d f x dj , (6) 

with x the structure parameter of the quadrupole op- 
erator of the IBM. For the Hamiltonian (JSJ) the energy 
surface in the large N limit reads (see for instance Eq. 
(6) in Ref. [2l| taking the large N limit and 7 = 0, this 
last condition is not a restriction since it is known that 
the IBM-1 including up to two-body interactions does 
not produce triaxiality) 

Eibm(N,(3) = (N,0\Hi BM \N,0) 

= Ar 7J^)l{5^-4-4y|/3x(^l) + 

P 2 [x+ 2 - X 2 ( X -l)]}- (7) 

The equivalence between (0} and (0) is readily established 
in the thermodynamical limit setting y = —y /l 2/7x- This 
correspondence can also be obtained from (jSJ if one is 
restricted to the spherical /x = component of the d- 
bosons. In that case 

Q = s^d + dls + X (2020|20)«4do, (8) 

with the Clebsch-Gordan coefficient (2020 1 20) = -y/2/7. 
It should be clear that this equivalence is only valid in 
the thermodynamical limit and restricted to the scalar 
excitations. Differences between both models are due 
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to the fluctuations induced by the 7-vibration over the 
mean field. The present scalar model does not include 
K = 2 excitations and consequently can not be used to 
study 7-excitations. But, apart from that, it mimics the 
usual representation of the IBM in a Casten triangle as 
shown in Fig. where any point in the triangle can be 
described in polar coordinates (p, 9). The radial variable 
p is related to x (basically p = 1 — x) and the angular 
variable 9 is related to y (9 — (1 — \/2y)-K/3) (or \ m the 
IBM). Thus, changing the control parameters (x,y) (or 
x, x m IBM) one can reach any point in the triangle and, 
consequently analyze any trajectory in it looking for crit- 
ical phenomena. In this respect, the present model allows 
to perform numerical calculations within the Casten tri- 
angle for large TV values, which are not posible with the 
usual IBM codes and are needed to check finite-size cor- 
rections of different observables as we will show in the 
next section. 
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FIG. 2: (Color online) Energy surfaces for the Hamiltonian 
(A in the large N limit QJ for y = l/y/2 (full line) and 
y — (broken line) and for different a;-values as a function of 
the deformation parameter f3. The limits of the coexistence 
region are shown in the panels x = 0.82036... (spinodal) and 
x — 4/5 (antispinodal) . The critical point is represented in 
the middle (left) panel (x = 9/11). The two panels on the 
right show two cases in the coexistence region. 

We would like to emphasize here that the st boson 
model has a natural order parameter, the expectation 
value of the i-boson number operator n t . This is not the 
case in the spin representation of the model • O n this 
regard, the Schwinger mapping could be used to translate 
the number operator n t into spin operators suggesting 
an order parameter for the latter model. In the coherent 
state representation (N — > 00) the expectation value of 
the i-boson number operator is simply given by (nt) /N = 
f3 2 / (1 + f3 2 ). For i=l, this order parameter vanishes 
whereas for x = it is given by 1 2 At 

6 J 4-H/ 2 + Vy 2 (4+2,2) 

the transition point x c , both minima lead to the same 
energy and the jump of the order parameter is y 2 /(4 + 
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FIG. 3: (Color online) Usual representation of the IBM phase 
diagram as a Casten triangle. The three dynamical symme- 
tries are located in the vertices. Critical points e(5) (second 
order) and x(5) (first order) with the coexistence region are 
also included. Variation of the control parameters, either in 
IBM or in the present model, allows to explore the whole 
model space represented by the triangle. 



y 2 ) which vanishes, as expected, for y = 0. We have 
displayed in Fig. 0]the behaviour of the order parameter 
as a function of x for two different values of y. 
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FIG. 4: (Color online) Variation of the order parameter as a 
function of x near the critical point. Full line are the analyt- 
ical results for N — 00 and circles are numerical results for 
N — 5000. The jump is clearly observed for y ^= 0. 



III. FINITE SIZE CORRECTIONS: 
MEAN FIELD 



BEYOND 



For y = 0, the st model corresponds, up to some re- 
definitions of the parameters, to the LMG model which 
has been recently shown to display nontrivial finite-size 
scaling properties at the transition point 0, 0, 0] . As 
explained above, the transition for y 7^ is of first order 
and it is thus natural to expect some qualitative changes 
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for the finite-size corrections. To analyze them, we follow 
the same route as those detailed in Ref. 0] for y = 0. 

For sake of clarity, we briefly sketch now this approach 
in the symmetric phase and refer to Ref. |17| for more 
details. For x > x c , we map the one-body operators 
using the Holstein-Primakoff boson expansion onto 
a single 6-boson 

th = tfb, (9) 
tU = N 1 / 2 tf(\-n b /N) 1 l 2 = {sh)\ (10) 
sU = N-n b . (11) 

with (nb)/N <C 1. At order p, in the symmetric phase, 
the Hamiltonian schematically reads: 

H =ii§i+ o ^i NP+i )- ( i2 ) 

i=0 

We make use of the canonical transformation method 
which allows to diagonalize exactly the Hamiltonian or- 
der by order. At order (1/N)°, H is simply a quadratic 
form and is thus simply diagonalized via a Bogoliubov 



transform. As explained in Ref. |l9j, diagonalizing H 
beyond this order only requires to solve a linear set of al- 
gebraic equations. More concretely, let us introduce the 
a-boson through: 

j=0 k,l 

where the are coefficients to be determined such that 
H expanded at order 1/N P is a polynomial function in n a . 
Further, one also has to impose the bosonic commutation 
rules which reads [b, tf] = 1. At each order, the content in 
and a is given by a simple power counting analysis. For 
p = 0, this transformation is nothing but the Bogoliubov 
transformation = o^a^ +a^\a and we thus have two 
coefficients to determine which are solutions of quadratic 
equations. At order p = 1, one must include linear and 
cubic terms and we have six coefficients ot^\ to find which 
are now solutions of a linear set of equations. One then 
gets, at order 1/N, the following diagonal form of the 
Hamiltonian: 



H = 



1 r 



x + E(x)^ 2 



-8 + 26x - 20x 2 + y 2 (-2 + 6x - 3x 2 ) 32 - 80x + 50x 2 + y 2 (8- Ux + 5x 2 ) 

h x 



x 2 (x- 1) 
2N 

nj-ix^ + ^pl 



S(x) 5 / 2 

lQx + y 2 (-2 + x) -16 + b2x - AOx 2 + y 2 (-4 + 16x - ICte 2 



S(x) 3 / 2 



E(x) 2 



: n 



_, _ x 2 {x - 1) -8 + 26x - 20x 2 + y 2 (-2 + 8x - 5x 2 ) 



N 



E(x) 2 



0(l/N 2 ), 



(14) 



with S(x) = x(5x — 4). As can be readily seen, the 
ground state energy as well as the gap are regular func- 
tions provided x > 4/5 which is the value of the critical 
point for y = 0. For nonvanishing y, it means that for 
x > x c > 4/5, there is no divergence in these quantities. 
One thus naively expect some 1/iV corrections from this 
side of the transition. 

The same approach can obviously be applied in the 
broken phase to diagonalize H at the same order. It 
then requires the exact knowledge of /?o which is solution 
of the cubic equation dE(N, (3)/d(3 = 0. If, for arbi- 
trary < x < x Cl it is a complicated expression of x 
and y, the nonvanishing solution for x — x c simply reads 
Po = y/2 and thus allows for simplifications. The main 
result of this study is that the Hamiltonian expanded 
at order 1/N is exactly the same at the critical point 
x = x c for Po = and /?o = y/2. Of course, at order 
N, one finds H = which simply confirms the existence 



of a first order quantum phase transition. The most in- 
teresting point is that the finite-size corrections are also 
the same for both values of /3o- Physically, it means that 
for x = xf the gap is A = E(x c ) 1/2 + 0{1/N) but for 
x = x c , the gap must vanish in the large N limit. In Fig. 
|S]we show the behaviour of the gap in the region around 
the critical point, x c = 9/11, for y = l/y/2. As can be 
seen gaps in the symmetric and the broken phases are 
indeed equal at x — xf. However, exactly at the critical 
point both states are degenerated giving rise to a zero 
gap as indicated by the black dot at x c = 9/11. In the 
inset, we show the behaviour of the gap at the critical 
point as a function of the boson number N, confirming 
the predicted exponential A ~ e~ AN . The exponent A is 
a function of y that vanishes in the limit y = where we 
expect A ~ iV -1 / 3 20]. This exponential decay clearly 
indicates that the 1/N expansions from the symmetric 
phase (/3q — 0) and from the broken one (/3q = y/2) are 
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FIG. 5: (Color online) Behaviour of the gap as a function of 
x near the critical point for y — 1/ vz. The full bold line and 
the full dot at x c — 9/11 correspond to the analytical results. 
Open dots are the numerical calculation for N — 5000. In the 
inset, the exponential decrease of the gap at the critical point 
is shown for several values of y. 
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FIG. 6: (Color online) Same as Fig. EJbut for y = 0. The full 
line is the analytical result for the gap, while the open dots 
correspond to the first two excited states for N = 5000. 



the same at all orders since any difference would imply a 
power-law decrease. As will be discussed in a forthcom- 
ing publication for the standard IBM model, we empha- 
size that this matching of the 1/N expansions is not a 
rule but rather an exception 0. In Fig. El the gap for 
y = is shown in the region around the critical point. 
The full line corresponds to the analytical expression and 
vanishes at the critical point. In this case, there are no 
parity breaking terms so results discussed in Ref. [l^ 



are recovered. For x > 4/5, the system is in the symmet- 
ric phase with positive parity. The gap is related to the 
energy of the first excited state corresponding to a one 
phonon state with negative parity. The second excited 
state is a two phonon state with positive parity. The ex- 
citation energies for these two states for N = 5000 are 
plotted in Fig. [5] as circles joined by a (green) line. For 
x < 4/5 the system is in the broken phase. In the limit 
x = there are two sets of states of different "parity" de- 
generated (± parity doublets) as displayed by the lowest 
(green) circle-line in Fig. [5] In this phase the nonvanish- 
ing gap is related to the excitation energy of a different 
phonon which corresponds to the upper (green) circle- 
line. It is equivalent to the f3 excitation in IBM (there 
are no 7 excitation in this scalar model). In the limit 
N — ► 00 both gaps (from the spherical and from the bro- 
ken sides) go to zero at the critical point, although they 
remain finite at finite N as clearly seen in Fig. 



IV. SUMMARY AND CONCLUSIONS 

To summarize, we have introduced a two-level model in 
terms of two scalar bosons s and t, whose phase diagram 
in the thermodynamical limit is the same as the IBM 
phase diagram when described by the QCH. We have 
exactly diagonalized the Hamiltonian in the symmetric 
phase up to order (1/7V) 1 and obtained the finite-size 
corrections for the ground state energy, the gap and the 
order parameter. These corrections were tested against 
numerical results for large systems. Though the phase 
diagram of this st model and the QCH coincide in the 
large N limit, beyond the mean- field level the QCH may 
have different properties due to the dynamics of fluctua- 
tions and correlations in the full (/?,7) space 0, |U l23l . 
In particular, while our st model is fully integrable |24| 
the IBM is chaotic for \ ^ j2f| with the exception of 
the ST/ (3) symmetry limit. Extension of the formalism 
presented in this paper to the IBM will be the subject of 
a forthcoming publication • 
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